Routes to gelation in a clay suspension 
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The gelation of water suspension of a synthetic clay (Laponite) has been studied by dynamic 
light scattering in a wide range of clay weight concentration {Cw = 0.003 0.031). At variance 
with previous determination, indicating a stable liquid phase for < ~ 0.015 ~ 0.018, we find 
that the gelation takes actually place in the whole examined Cw range. More importantly, we find 
that marks the transition between two different routes to gelation. We hypothesize that at low 
concentration Laponite suspension behaves as an attractive colloid and that the slowing down of the 
dynamics is attained by the formation of larger and larger clusters while at high concentration the 
basic units of the arrested phase could be the Debye Huckel spheres associated to single Laponite 
plates. 

PACS numbers: 64.70.Pf, 82.70.Dd, 78.35.+C, 61.20.Lc 
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Colloidal systems are ideal benchmarks for studying, 
by optical microscopy and light scattering, equilibrium 
slow dynamical processes and the formation of non er- 
godic arrested state of matter. In this context sterically 
stabilized colloidal systems have been often used as model 
for hard spheres, to study the fluid-crystal transition and 
the glass formation^. For suspensions of hard spheres 
the latter transition occurs at relatively high packing 
fraction {(p > 0.58). When screened charged interactions 
are present, crystal and disorder arrested phases can be 
formed also at very low volume fraction ((p « 0.01). Re- 
cent works 1^ y, y, y| focusing on short-range attractive 
colloids have attempted a connection between the gel and 
the glass arrested state of matter. In these systems, a re- 
entrant glass line, two kind of glasses (attractive and re- 
pulsive) and a glass-glass line have been predicted 
and experimentally observed Correlations between 

the dynamical behavior of gels and glasses suggest that 
a common understanding of these two disordered forms 
of matter may emerge. In the case further complicated 
where short-ranged attractive interactions are comple- 
mented by weak repulsive electrostatic interactions, the 
gel formation process can be fully modeled as a glass tran- 
sition phenomenon 7^. In the last case, the gel phase is 
stabilized by the competition of the short range attrac- 
tion and the long range repulsion. Despite these recent 
progresses, a deeper comprehension of the still puzzling 
liquid-gel transition in colloidal system is requested. 

Aqueous suspensions of charged platelike colloids have 
been the subject of intense experimental and theoretical 
investigations. The synthetic Hectorite clay Laponite is 
a perfect model for these charged platelike colloids. It 
is in fact composed of nearly monodisperse, rigid, disc- 
shaped platelets with a well defined thickness of 1 nm, an 
average diameter of about 30 nm and a negative surface 
charge of a few hundred e. In spite of intensive research 
on Laponite suspensions motivated by its important in- 
dustrial applications, there is no general agreement about 



the mechanism of gelation; both attractive and repulsive 
forces are claimed to be responsible for the gel formation 
attributed by different authors to Wi ener glass transition 
||8||] , frustrated nematic transitio n [p. lloL [ill , micro- 
segregation gelation |l5l ll6L ll^ etc. Hence, 
studies of the aggregation in Laponite suspensions will 
give an important contribution to the understanding of 
the gel formation process in systems where repulsive long 
range and attractive short range interactions compete. 

According to the phase diagram obtained from 
Mourchid et al. at ionic strength below / = 10~^ 
M Laponite suspensions can be in two different physical 
states. Low concentration suspensions [Cw < C^{r)) 
form stable, equilibrium fluid phase. Higher concen- 
tration suspensions (C^, > C*,(/)) are initially fluids 
but experience aging and pass into a gel phase after a 
time that depends on the clay amount. Several dynam- 
ical scattering studies have investigated the aging pro- 
cess of samples belonging to this high phase region 
[^JJj JJ., _20j 21], while only two recent papers ^,,21] re- 
port about the lower concentration region. While the first 
paper fl^ reports about the formation of a gel that does 
not flow when tilted at low laponite concentration {C'w—1 
%, 1= 5 mM NaCl), the authors of the second paper |2]| 
find that samples at low concentrations (Cu.=(0.1-^1) %, 
/ < 5 • 10~*M)are stable liquid suspensions. For this 
reason and to try to understand the gelation mechanism 
for this charged colloidal system we have performed ac- 
curate and systematic measurements on several samples 
in the low concentration region and at higher clay con- 
centrations. 

In this letter we present a dynamic light scattering 
study at increasing concentrations from 0.3 to 3.1 wt % 
at ionic strength between / ~ 10"^ M and / ~ 10"'^ 
M where C* ^ 1.5 ^ 1.8% g. A new gel region, at 
Cw < C^,) has been identified and the gelation processes 
in the two gel phases (at Cw < and at C„ > C^) have 
been accurately studied. Two different routes to gelation 
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for the two regions have been found and a possible expla- 
nation for the processes involved in the building of the 
arrested phases has been proposed. 

Laponite is a hectorite synthetic clay manufactured by 
Laporte Ltd, who kindly supplied us the material. Par- 
ticular attention has been devoted to samples preparation 
to avoid dissolution of Laponite platelets that occurs if 
the samples are exposed to air contamination due to the 
presence of atmospheric CO2 113, 113 and can seriously 
affect the measurements. For this reason the whole pro- 
cedure has been performed in a glove box under N2 flux 
and the samples have been always kept in safe atmo- 
sphere during and after sample preparation. The powder 
was firstly dried in an oven at T=400 K for 4 hours (up 
to 20% of the as received powder weight is due to ad- 
sorbed water), then the powder is dispersed in deionized 
(pH=7) water, stirred vigorously until the suspensions 
were cleared and then filtered through 0.45 /im pore size 
Millipore filters. For same of the samples the pH has been 
measured by a Crison Glp 22 pHmeter, the final value, 
reached after 10-^20 min from the laponite dispersion in 
water, was in the range pH=9.8-;-10.0. The final value 
of pH=10, neglecting the contribution of Na counterions 
released from some of the platelets to the solution, would 
fix the ionic strength around 10~''M. Assuming that 50 
counterions are released by each laponite platelet fisf . 
the ionic strength never exceeds ~ 10~^M for all the in- 
vestigated samples. The fact that laponite dissolution 
does not take place in our samples is confirmed by NMR 
spectroscopy [2^ that checked continuously the amount 
of Na and Mg ions present in the solutions. These mea- 
surements show in fact that, both in the low and in the 
high concentration ranges, the fraction of solvated Na 
ions does not change with time, and that the amount of 
solvated Mg ions are below the instrumental sensibility. 
These observations indicate that the chemical reaction of 
Laponite dissolution [UH^ is not affecting our samples. 
Futhermore, in none of the samples studied we observed 
any evidence of sedimentation also several months after 
their preparation as has been checked by direct sample 
inspection and by performing test measurements of the 
correlation function at different heights in the cells. The 
starting aging time (tu)=0) is defined as the time when 
the suspension is filtered. This sample preparation pro- 
cedure is similar, but not identical, to that already used 
in previous works 0, 0, HO, 0| . Some authors do not 
report information on the drying of the powder, others 
allow the suspension to get in contact with atmosphere, 
and different details can lead to small differences in the 
sample behavior on aging. In addition, the determina- 
tion of the sample concentration is also problematic, as 
the filtering procedure can alter the actual concentration 
if large clusters are present in the suspension after the 
stirring phase. To minimize this effect, we have prepared 
a large amount of concentrated solution, then -after the 
filtering- we have diluted this parent sample at the de- 



sired concentrations. In this way, in the present mea- 
surements, the relative concentrations are well defined, 
even if the absolute value in Cyj can be affected by a sys- 
tematic error, that we estimate in the range AC„,/Ctu ~ 
0.1. 

Dynamic light scattering measurements were per- 
formed using an ALV-5000 logarithmic correlator in com- 
bination with a standard optical set-up based on a He-Ne 
(A = 632.8 nm) 10 mW laser and a photomultiplier de- 
tector. The intensity correlation function was directly 
obtained as 52(9,*) =< Iiq,t)I{q,0) > / < /(g,0) >^, 
where q is the modulus of the scattering wave vector 
defined as q = {Ann/ X)sin{d /2) {9—90° in the present 
experiment). For the low concentration samples (C^ < 
1.5%) the measurements have been repeated once a week 
for a long period of time (up to four months for the 
lowest concentration), while for Cw > 1.5% the sample 
was left in the scattering system and photocorrelation 
spectra were continuously acquired until the gelation was 
reached. 
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FIG. 1: Evolution of the measured intensity correlation func- 
tions (symbols) and corresponding fits with Eq. Q (continu- 
ous lines) for two different Laponite suspensions at the indi- 
cate concentrations at different waiting times t^,. The curves 
are measured at increasing waiting times that (from left to 
right) are t™=288, 792, 1128, 1464, 1632, 1800, 1968, 2136, 
2328, 2640 hours for sample (A) and t„=6, 30, 54, 78, 102, 
126, 150 hours for sample (B)(for which only data in the er- 
godic phase are reported). 
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As an example, correlation functions at increasing ag- 
ing times tw for two different samples at low- (Fig. lA 
- Cm = 1.1 wt % < ) and high- concentration (Fig. 
IB Cu, = 2.5 wt % > C^), are reported in Figure 1. As 
it is evident from the figure both the samples are doing 
aging, the dynamics is in fact becoming slow and slow for 
increasing waiting time t^,- This behavior was expected 
for the higher concentration sample of Fig. IB, several 
measurements 0, 0, 0, HO, 0| , indeed, have shown the 
existence of aging and a sol-gel transition for long enough 
waiting times. The aging process was instead unexpected 
for the lower concentration sample of Fig. lA where the 
liquid state is predicted from the phase diagram , and 
confirmed by recent dynamic light scattering measure- 
ments [i^ where, perhaps, the authors did not wait for 
a time long enough to observe significant aging. From 
Fig. lA it is instead evident that the system is aging and 
for the longest waiting time reported (ttu=2640 hours) 
there is a qualitative change in the correlation function: 
a crossover between a complete and an incomplete de- 
cay. This behavior is the indication of a strong ergodic- 
ity breaking, signature of a sol-gel transition, as already 
observed in (3j uM fo^ samples in the higher concentra- 
tion region and in for a sample at = 1% with 
I — 5 ■ 10~^M . Our dynamic light scattering measure- 
ments show that for all the samples studied, down to the 
lowest concentrated one (C^ = 0.3%), there is a typical 
waiting time, , increasing with decreasing clay concen- 
tration Cw, at which the system undergoes the gelation. 
Depending from the initial concentration this gelation 
time can vary from hours to several months. 

Figure 1 shows also that the correlation functions de- 
cay following a two steps behavior, i. e. there are two 
different relaxation processes, a fast and a slow ones. For 
this reason the fitting expression should contain two con- 
tributions. In this case the squared sum of an exponential 
and a stretched exponential function is used, as already 
reported in Ref.[20|: 

g2iq,t) -1 = 6 (ae(-*/"i) + (1 - a)e(-(*/"^)''))' (1) 

where b represents the coherence factor. The stretched 
exponential, instead than another simple exponential 
term, is used since it has been found to give a good de- 
scription of the slow relaxation process in glassy systems. 
The fitting expression well describe the photocorrelation 
spectra for all the aging times in the liquid (ergodic) 
phase and for all the investigated concentrations with 
the relaxation time ti associated to the fast dynamics 
and the relaxation time T2 and the stretching parameter 
f3 that describe the slow part of the autocorrelation func- 
tion. The fits are shown as full lines in Fig. 1 and the 
resulting chi square is always within its standard devia- 
tion. 

Another important peculiarity of the correlation func- 



tions, that can be directly observed in the raw data re- 
ported in Fig. 1, is that the aging process evolves dif- 
ferently for the lower and higher concentration samples. 
While in fact the initial decay of the correlation functions 
-fast dynamics, characterized by the correlation time ri- 
remains constant with waiting time for the higher con- 
centration sample (Fig. IB) it seems to become slower 
with increasing waiting time for the lower concentration 
sample (Fig. lA). This different behavior is confirmed 
from the results of the fits, that also show another major 
difference between the samples at low and high concen- 
trations. 

In the following we will report the tw and Cw depen- 
dence of the fitting parameters, in particular we focus on 
the behavior of the parameters of the slow decay: the 
relaxation time T2, the stretching exponent /3 and the 
"mean" relaxation time r^: 

T™ = T2ir(i). (2) 



where F is the usual Euler gamma function. 
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FIG. 2: Waiting time dependence of Tm (see Eq. I^J, i. e. of 
the average relaxation time of the slow decay of the correla- 
tion functions as those reported in Fig.0 As an example, the 
dependence of Tm is reported for the three indicated con- 
centrations. Continuous lines are fits of the data with Eq.|3 

As an example, the values of r„i obtained for three dif- 
ferent concentrations are reported in Fig. 2. It is evident 
the common behavior of Tm that seems to diverge at a 
given tu, = i.e. when the gelation occurs. This clearly 
shows that, as already seen from the raw correlation func- 
tions, also the samples at concentrations lower/equal to 
1.5 wt % are actually doing aging and undergo a gelation 
transition. To have more information about this gelation 
process we represent the aging time (t^,) dependence of 
mean relaxation time (r™) with the law: 
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FIG. 3; Concentration dependence of the divergence time t'ij^ 
(A) and of the B parameter (B), entering in Eq. |3 The 
dashed lines are guides to the eyes. The shadow area indicates 
the region where the liquid-gel transition was supposed to be 
according to Ref. Q 



-r™ = To exp [b J"^ ) (3) 

We do not have any rationale behind such a law, it can 
be considered a generalization to long waiting time of the 
exponential growth with tw observed for T2 by Abou et 
al. |20l | in the high concentration samples. Equation |21 
-reported in Fig. 2 as full lines- well describe the mea- 
surements. The most significative parameters of the fits, 
and B, for all the studied concentrations are shown in 
Fig. 3. Here the vertical dotted region indicates the range 
of concentrations that, in the phase diagram 0, would 
mark the transition from the "liquid" to the gel phase 
in the range of / ~ 10"^ - IQ-^'M. The results of the 
fit indicate that -that can be considered as the time 
at which the gelation actually occurs- is continuously de- 
creasing with increasing clay concentration, without any 
evident discontinuity in correspondence of the "transi- 
tion" region. The concentration dependence of the pa- 
rameter B, which measures how fast t,„ approaches the 
divergence, is shown in Fig. 3B. This parameter is almost 
constant for all the samples in the low concentration re- 
gion, while shows a large discontinuity on passing in the 
higher concentration region. It is important to note that 



this jump takes place in a region that encompass the sup- 
posed "liquid-gel transition" region. 




FIG. 4: The waiting time dependence of the Tm (A) and /3 
(B) are reported as a function of the scaled variable t^/t^ for 
some of the investigated concentrations: ■ = 0.3%, • = 0.5%, 
A = 0.6%, ▼ = 0.8%, ♦ = 0.9%, + = 0.9%, X = 1.0%, 
* = 1.1%, □ = 1.2%, o = 1.3%, A = 1.4%, V = 1.5%, 
□ = 2.2%, El = 2.5%, B = 2.8%. The data for C„ < 1.5% 
collapse on a single master curve. 

The fact that the value of the B parameter is almost 
constant for all the low concentration samples is an in- 
dication of the existence of a scaling law. Indeed, if we 
plot the quantities Tm{tw) and P{tw) as a function of 
tw/tw ^-ll the different Cw data with almost constant B 
should collapse on a single master curve. The Tm and 
the (3 parameters in function of the normalized waiting 
time are reported in Fig. 4A and 4B respectively. As 
expected, all the data for both the Tm and the (3 param- 
eters of lower concentrations collapse, within their sta- 
tistical uncertainties, on a single curve while the data at 
higher clay concentrations have a different behavior. We 
already observed from the direct comparison of Fig. lA 
with Fig. IB that the aging process is qualitatively differ- 
ent for the low and high concentration samples. Figures 3 
and 4 quantify this difference in the physical properties 
characterizing the aging phenomenon in the two different 
concentration regions. 

In conclusion the present observations indicate that 
the stable phase of Laponite suspensions in pure water 
at Cw ^ C*^ -expected to be liquid according to previ- 
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ous studies [91- a ctually is an arrested phase. Recently 
Nicolai et al. [13 proposed a reyisitation of the phase di- 
agram and suggested that the equilibrium state in the low 
concentration region is not liquid but a very fragile gel 
that takes a long time to appear 0| . Probably this " long 
time" is the reason why previous measurements indicated 
this phase as liquid. Our measurements clearly show not 
only that at all the investigated concentrations the stable 
phase is an arrested phase but also that the aggregation 
process is basically different for the samples at < 
and those at ^ C'^ . This is evident from the behavior 
of Tm and (5 reported in Fig. 4 and from the discontinuity 
of the B parameter across the C^, = region as shown 
in Fig. 3. In this sense the "liquid-gel" transition line of 
the phase diagram is not a real liquid-gel transition 
but rather seems to indicate a sort of "fragile gel-gel" 
transition. 

The origin of the different routes towards an arrested 
phase observed for low- and high-concentration Laponite 
suspensions calls for an explanation that goes beyond the 
aim of the present work. At the level of speculation, it is 
worth to recall that the microscopic interaction between 
Laponite plates is due to a screened Coulomb interac- 
tion, which can be modeled by a Yukawa-like repulsion 
at long distances and a quadrupolar electric and/or a van 
der Waals terms at short distances [i^ ■ This com- 
petition of short-range attraction and long range repul- 
sion resembles that recently proposed to describe the 
phenomenology of colloidal gels. In Ref. it has been 
suggested that the gelation taking place in attractive col- 
loidal suspensions at very low concentration involves the 
growth of larger and larger clusters (driven by the short 
range attraction). The glass transition of these clusters, 
driven by the long range repulsion, is responsible for the 
arrest of the dynamics and the formation of a gel phase. 
The authors prove that the gel is essentially a Wigner 
glass, composed by clusters of colloidal particles. We can 
speculate that low concentrations Laponite suspensions 
behave as the attractive colloids in Ref. • This scenario 
is in agreement with the fact that the long time needed 
to form the arrested phase in the low concentration range 
is spent by the system to build up the clusters. Also at 
high Laponite concentration the arrested state could be 
a Wigner glass but in this case the packing fraction of the 



Debye Huckel sphere associated to each platelets reaches 
values as high as 0.43 and the single laponite platelet 
would be the elementary constituent of the Wigner glass. 
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